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Abstract 

We consider J\f = (1, 1) super Yang-Mills theory in 1+1 dimensions with 
fundamentals at large- N c . A Chern-Simons term is included to give mass to the 
adjoint partons. Using the spectrum of the theory, we calculate thermodynamic 
properties of the system as a function of the temperature and the Yang-Mills 
coupling. In the large-iV" c limit there are two non-communicating sectors, the 
glueball sector, which we presented previously, and the meson- like sector that we 
present here. We find that the meson-like sector dominates the thermodynamics. 
Like the glueball sector, the meson sector has a Hagedorn temperature Th, and 
we show that the Hagedorn temperature grows with the coupling. We calculate 
the temperature and coupling dependence of the free energy for temperatures 
below Tjj. As expected, the free energy for weak coupling and low temperature 
grows quadratically with the temperature. Also the ratio of the free energies at 
strong coupling compared to weak coupling, r s _ w , for low temperatures grows 
quadratically with T. In addition, our data suggest that r s ^ w tends to zero in the 
continuum limit at low temperatures. 



1 Introduction 



The free energy of M = 4 super Yang-Mills (SYM) theory at a large number of colors 
N c is larger at strong coupling by a factor 4/3 compared to weak coupling [H |2j. The 
weak-coupling result is calculable in perturbation theory, while the strong-coupling re- 
sult can be derived from black-hole thermodynamics. In the light of this finding, one 
can ask if other SYM theories exhibit a similar behavior. An analytic calculation of 
both the strong and the weak coupling limit of a field theory is generally not possi- 
ble, although there have been a number of proposals for methods to obtain solutions 
of finite-temperature supersymmetric quantum field theories [3]. We will use a nu- 
merical approach which is based on Supersymmetric Discrete Light-Cone Quantization 
(SDLCQ) [HE], thereby preserving supersymmetry exactly. Currently, this is the only 
method available for numerically solving strongly coupled SYM theories. Conventional 
lattice methods have difficulty with supersymmetric theories because of the asymmetric 
way that fermions and bosons are treated, and progress [6] in supersymmetric lattice 
gauge theory has been relatively slow. 

Previously we calculated the thermodynamic properties of pure glue M = (1,1) 
SYM theory in 1 + 1 dimensions [7j. Here we extend the calculations to include a sector 
with fundamental partons. In the large- N c limit the bound states in this sector of the 
theory are chains in color space with a fundamental parton at each end. The links in 
the chain are adjoint partons. Bound states of this type will be called mesonic because 
they have two fundamental partons, whereas solutions with only adjoint partons will 
be called glueballs. We have also extended the calculation to include a Chern-Simons 
(CS) term, which gives mass to the adjoint partons. 

The mesons and glueballs constitute two sectors of the same theory; both contribute 
to the thermodynamics. In the large-iV c limit the sectors decouple. Due to the cyclic 
redundancy of single-trace glueball states, there are many more meson states, which 
are in turn likely to dominate the thermodynamic properties of the system. In pre- 
vious work on the glueball sector [7] we found that the system possesses a Hagedorn 
temperature. 

Recall that SDLCQ makes use of light-cone coordinates, with x + = (x° + x 3 ) / y/2 
the time variable and p~ = {p° — p 3 )/\/2 the energy. One must be careful in defining 
thermodynamic quantities on the light-cone. It seems natural to define the partition 
function [8] on the light-cone as e~^cP . However, as shown by Alves and Das [9] , the 
above prescription leads to singular results for well known quantities which are finite 
in the equal-time approach. Their argument is based on the fact that using e~^cP as 
the partition function implies that the physical system is in contact with a heat bath 
that has been boosted to the light-cone frame. However, this is not equivalent to the 
physics of a system in contact with a heat bath at rest. This can be realized in a more 
direct way by noting that, since the light-cone momentum 

(p° + p 3 ) 

p ^—fT 

is conserved, the partition function must include the conserved quantity and is of the 
form 
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where /i is the chemical potential corresponding to the conserved quantity. The phys- 
ical interpretation of the chemical potential is that of a rotation of the quantization 
axis. Thus \x = 1 corresponds to quantization in an equal-time frame, where the heat 
bath is at rest and the inverse temperature is j3 — a/2/3 lc , and /i ^ 1 corresponds to 
quantization in a boosted frame where the heat bath is not at rest. Thus \x corresponds 
to a continuous rotation of the axis of quantization, and /i = would correspond to 
rotation all the way to the light-cone frame. 

A rotation from an equal-time frame to the light-cone frame is not a Lorentz trans- 
formation. It is known that such a transformation can give rise to singular results for 
physical quantities. This appears to be consistent with the results found in [9]. A 
number of related issues have been extensively discussed by Weldon [10]. The method 
has also recently been applied to the Nambu-Jona-Lasinio model 

The difficulties are avoided if we compute the equal-time partition function Z = 
e~ l3p ° , as was proposed much earlier by Elser and Kalloniatis [12]. The computation 
may, of course, still use light-cone coordinates. Elser and Kalloniatis did this with 
ordinary DLCQ [131 EH] as a numerical approximation to (l+l)-dimensional quantum 
electrodynamics. 

Here we will follow a similar approach using SDLCQ to calculate the spectrum of 
M = (1,1) super Yang-Mills theory in 1 + 1 dimensions [T5]. Though this calculation is 
done in 1+1 dimensions, it is known that SDLCQ can be extended in a straightforward 
manner to higher dimensions [161 El HE]- As is customary, we will assume that the 
single-trace bound states of our large- N c approximation are single-particle states. 

We have discussed the SDLCQ numerical method in a number of other places, and 
we will not present a detailed discussion of the method here; for a review, see [5J. 
For those familiar with DLCQ [131 EH], it suffices to say that SDLCQ is similar; both 
impose periodic boundary conditions on a light-cone box x~ G [-L, L] and have discrete 
momenta and cutoffs in momentum space. In 1 + 1 dimensions the discretization is 
specified by a single integer K = (L/n)P + , the resolution [T3], such that longitudinal 
momentum fractions are integer multiples of 1/K. However, SDLCQ is formulated 
in such a way that the theory is also exactly super symmetric. Exact supersymmetry 
brings a number of very important numerical advantages to the method; in particular, 
theories with enough supersymmetry are finite. We have also seen greatly improved 
numerical convergence in this approach. 

The calculation of thermodynamic quantities requires summing over the spectrum 
of available states, which we represent by a density of states (DoS). We will use a 
new numerical approach to estimate the density of states. The new approach is more 
efficient than the method used in previous work [7], because it allows us to extract 
the density of states without fully diagonalizing the Hamiltonian, a computationally 
challenging task. This innovation will enable us to pursue calculations at higher values 
of the resolution K. 

We find that the two-dimensional SYM theory with fundamentals and a Chern- 
Simons term exhibits a Hagedorn temperature Th [19], and calculate Th for several 
values of the resolution K and Yang-Mills coupling g. Extrapolating to the continuum 
limit, we obtain Th as a function of the coupling. The Hagedorn temperature is used 
as an upper limit for the temperatures we can use to calculate the thermodynamic 
properties of the system. 

In Sec. [2] we provide a review of the formulation of super Yang-Mills theory with 



2 



fundamental matter and a Chern-Simons term in 1+1 dimensions. In Sec. [3] we dis- 
cuss some of the properties of the SDLCQ spectra in some limiting cases and provide 
comparisons between glueball and mesonic sectors of the theory. The discussion in 
Sec. H] presents the methods for estimating the density of states and the Hagedorn 
temperature. In Sec. [5j we summarize our formulation of the thermodynamics and the 
formulae we use to calculate the free energy. We then present the numerical results for 
the free energy, which we obtained using the DoS approximation to the spectrum, at 
various values of SYM coupling g up to the Hagedorn temperature. Finally, in Sec. [6] 
we conclude by summarizing our results and the prospects for future work using these 
methods. 

2 Super Yang— Mills theory with fundamental mat- 
ter and Chern— Simons term 

2.1 Formulation of the theory and its supercharges 

We start by considering M = 1 supersymmetric gauge theory in 2 + 1 dimensions 
coupled to fundamental matter and a Chern-Simons three-form. The action is 

S2+1 = Sym + S{ 

.matter + S cs , (2.1) 

with 

Sym = j dx 3 Tr(~F^ + ^AT"D M A), (2.2a) 

^.matter = j dx* (D^D^ + MD^V - g [tf A£ + £ + Att] ) , (2.2b) 

Scs = Jdx 3 ^ (<T X {A,d u A x + jgApA v A x ) + 2AA^) . (2.2c) 

The SYM part of the action describes a system of gauge bosons A^ and their super- 
partners, the Majorana fermions A. Both fields are (N c x N c ) matrices transforming 
under the adjoint representation of SU(N C ); hereafter, unless indicated otherwise, we 
treat these fields as matrices, and thus we suppress the color indices k). Addition- 
ally, we have two complex fields, a scalar £, and a Dirac fermion all transforming 
according to the fundamental representation of the gauge group. In matrix notation 
the covariant derivatives and the gauge field strength are defined as follows: 

D fJ A = d tl A + ig[A fl ,A], = d^ + igA^, = 3^ +igA l &, 

(2.3) 

= d£ - iggAp, = 8^ - ig^Ap, = d^A v] + ig[A^ A,}. 
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The action (12 .ip is invariant under supersymmetry transformations parameterized by 
a constant two-component Majorana spinor e = (ei, £2) T ', £ = £ T T°: 

SA^ = ifr M A, 5A = i^r^e, 
5£=^>, ^ = ~*e, (2.4) 

where T^, the spinor generator of the Lorentz group, is written as 

V n» = = ie ^Y x ■ ( e - +2 = i). 

Using standard Noether techniques, we construct the spinor supercurrent corre- 
sponding to the above supersymmetric field variations 

eq»= -fr^I^Tr (AF a/3 ) + )-D*gEV + ^eT^D^ - ^leD^ + ^DjlY^ei. (2.5) 

For the remainder of the paper we assume that the fields are independent of the space- 
like dimension x 2 , i.e. <92(...)=0, thereby dimensionally reducing the theory to two 
dimensions. Thus the M = 1 supersymmetry in 2 + 1 dimensions is naturally expressed 
in terms of M = (1,1) supersymmetry in 1 + 1 dimensions. 

We will implement light-cone quantization, which means that initial conditions as 
well as canonical (anti-) commutation relations will be imposed on the light-like surface 
x + = const. In particular, we construct the supercharge by integrating the supercurrent 
(12. 5p over the light-like surface 

eQ = j dx~ Qer^r+Tr (AF aj3 ) + X -D_t} etf + ^eT +v D v * 

v (2.6) 

Note that because we have taken the fields to be independent of x 2 , the integration 
over this coordinate resulted in a constant factor, which rescaled our original fields. 

By choosing the following imaginary (Majorana) representation for the Dirac ma- 
trices in three dimensions: 

r° = a 2 , r^iai, r 2 = kr 3 , (2.7) 

the Majorana spinor field A is manifestly real, i.e. A* = A T . At this point it is 
convenient to introduce the component form for the spinors: 

A=(A,A) T , *=(^) T , Q=(Q+,Q-) T (2.8) 

In terms of this decomposition, the superalgebra is realized explicitly in its M = (1,1) 
form, namely 

{Q + ,Q + } = 2V2P + , {Q-,Q-} = 2V2P-, {Q+,Q-} = 0, (2.9) 
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where Q + (Q~) are left (right) Majorana-Weyl spinors, each characterizing the smallest 
spinor representation in 1 + 1 dimensions. 

To readily eliminate the nondynamical fields, we impose the light-cone gauge (A + = 
A_ = 0). In this case the supercharges can be read off (12. 6p and are given by 

Q + = \f dx~ ^Tr(Xd^A 2 ) + ^cL£V - ^t(9_£ - ^d-ij) + ^d-i/J%\ , (2.10) 

Q- = -L J dx~ (lr(\d_A-) - tfDrf + iD^i - ^=M^£ ~ £^)) • (2.11) 

Notice that the right-movers (ip) appear in the supercharge Q~ only in the total deriva- 
tive term. This is a consequence of the light-cone formulation, which singles out the 
non-dynamical fermion degrees of freedom, leaving in the expression only the physical 
spinor fields (A and ip). Among the equations of motion that follow from the action 
(12. ip . in the light-cone gauge, three serve as constraints rather than as dynamical equa- 
tions. Namely, for A and ip, respectively, we have 

dS\ = -^={[A\X}+ i^t _ m \ _ iftA) , (2.12) 

fl_^ = -^4V + ^. (2.13) 
While for A~ we obtain 

d 2 _A~ = gj, (2.14) 

with 

J = -i[A\ d^A 2 ] + -i={A, A} - i(0_Of f + i£(d-£ f ) + + -d^A 2 . (2.15) 

V2 9 

Note that the field A~ has to be eliminated from the supercharges, in favor of the 
physical degrees of freedom. This can be done by inverting (12. 14p . 

The only contribution from the Chern-Simons term enters into the supercharges 
via equation ( 12. 15ft . because 5Ccs oc <9 M (. . .) under the supersymmetry transformations 
(j2.4P ; the surface term is not shown here. The inclusion of a Chern-Simons term 
in our theory is important, since it effectively generates mass for the adjoint partons 
proportional to the coupling k. 

2.2 Bound-state eigenvalue problem 

The bound-state spectrum is obtained by solving the following mass eigenvalue equa- 
tion: 

2P + p-\<p) = V2P + (Q~) 2 \<p) 

= V2P+(g(Q SYM + Q- mattcI )+iKQ cs ) 2 \v) = M 2 |^>, (2.16) 
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where the various pieces of Q , after dropping the surface terms and eliminating the 
non-dynamical fields using the constraint (12.141) . may be expressed as follows: 



QSYM = ^j dX ~ (V' d ~ A ^ + 7f ^ ( 2 - 17a ) 

+ ^A 2 4> + ^A 2 ^, (2.17b) 

Qcs = ^=J dx-{d_A 2 )±-\ (2.17c) 

where a trace over color space is understood. 

The strategy for solving equation (12. 16[) is to cast it as a matrix eigenvalue prob- 
lem. This is achieved by employing a discrete basis where the longitudinal light- 
cone momentum P + is diagonal. The discrete basis is introduced by first discretizing 
the supercharg^j] Q~ and then constructing P~ from the square of the supercharge: 
P~ = (Q~) 2 / \^2. The two dimensional theory is compactified on a light-like circle 
(— L < x~ < L), and periodic boundary conditions are imposed on all dynamical 
degrees of freedom. This leads to the following field mode expansions: 

4(0, x -) = —^2— (a ld {n)e—-l L + o^n)^"^) , (2.18) 

* n=l 
1 oo 

6(0, a;") = —j= y~] —j= (ci(n)e~ innx ~ /L + c J t (n)e imnc ~ /i ) , (2.20) 

-. oo 

^(0, z") = ; V (di(n)e- imrx ~/ L + d\(n)e innx ~ /L ) . (2.21) 



2iV2L n=1 

In the above expressions^] we introduced the discrete longitudinal momenta k + = k 
as fractions nP + /K = nir/L; (n = 1,2,3,...) of the total longitudinal momentum 
P + , where K is the integer that determines the resolution of the discretization. The 
color indices were made explicit as well. Because light-cone longitudinal momenta are 
always positive, K and each n are positive integers. The number of constituents is thus 
bounded by K. The continuum limit is reached by letting K — > oo. 

The time direction in the light-cone formalism is taken to be the x + direction. 
Thus the (anti-)commutation relations between fields and their conjugate momenta are 
assumed on the surface x + = 0. Quantization is achieved by imposing the following 



■"■Note the relative phase between Qg Y M anc ^ Qcs- Qsym * s defined as Hermitian and Qq S is 
defined to be anti-Hermitian such that Q~ remains Hermitian. 

2 Thc inclusion of zero modes is beyond the scope of the present paper. 
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relations: 

[j§(0,x-),d-AUV,V~j\ = 1 (t^ik ~ 5{X ~ ~ V ~ ] ' (2 ' 22) 

{Aij(0, x~), A«(0, y")} = \/2 f <f«<f ifc - -^MwJ <S(aT - y"), (2.23) 

^(0,x-),9_0(0,r)] =i5 ij 5(x--y-), (2.24) 

{^(0,x-),^(0,y-)} = >/2£ y £(*~ -y~). (2.25) 

The above (anti-) commutators can also be expressed, with the help of equations 
(12.181) — (I2.2ip . in terms of creation-annihilation operators 



a i:j (k), altik') = \ bij{k), \ = ( 6^6 fl - -^6^ ) 5(k - k') 



(2.26) 



d(k), c)(k')\ = ]^{k),c){h!)\ = [di(k), d](k')j = [di(k), d){k')^=5 %J 5{k - k'). (2.27) 

The expansion of the supercharge Q~ in terms of creation and annihilation operators 
is a straightforward exercise. For instance, the decomposition of Qq S in terms of Fourier 
modes gives the following expression: 



Q 



cs 



2 5 / 4 v^F 



' n=l * ^ 



n)a ij {n / 



'2.28) 



Similarly, for the supercharge, Q f matter , that controls the behavior of the fundamental 
matter fields, we obtain 



Q { 



matter 



ig^/L\ ^ (n 2 + n 3 ) /_ t , 



2 5/4 



+ -^-^4( ra 2)4i( n l) rf i( n 3) + d\ (713)^(712)6^(711) + 6^(711)4(712)^(713) 



+ 4(^3)6^(711)^(712) 



c](n 3 )ay(n2)dj(ni) 



+ 0^(712)4(711)^(713) + 4(711)0^(712)^(713) + cj( n 3)4'( ra l) a ii( ra 2) 

+ ^-^==^a] i (n 2 )cl(n 1 )rf j (n 3 ) + dt (713)0^(712)^(711) 



+ 4( ri 3)ci(7ii)ay(7i 2 ) + c\(n 1 )a\ j (n 2 )d j (nz) j <5„ 3 ,„ 1+n2 



(2.29) 



Our computer code carries out these expansions automatically. 

Apart from supersymmetry, the theory we set out to explore possesses another 
symmetryjl which may be used to reduce the size of the Hamiltonian matrix we need 



3 Note that when n = 0, parity (V) is also conserved. 
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to produce and diagonalize. Namely, we have a Z2 symmetry T that is associated with 
the orientation of the large- N c string of partons in a state [201 f2Tj . It gives a sign when 
the gauge group indices are permuted 

a,ij(k) — >■ -a,ji(k), bij(k) — >■ -bji(k). (2.30) 

In this paper we will discuss numerical results obtained in the large- N c limit, i.e. 
terms of order 1/N C in the above expressions are dropped. Note that corrections on the 
order of 1/N C are expected to lead to interesting effects [22]; however, they are beyond 
the scope of this work. 

3 Meson and glueball spectra 
3.1 Limiting cases 

We first investigate the strong-coupling (g ^> k) and weak-coupling (g — > 0) limits 
of the theory. In the strong-coupling limit, we previously found [23] that there are 
approximate BPS glueball bound states with masses (squared) 

M 2 aB p S = n 2 K\ n = 2,3,...,(K-l). 

In the meson sector under investigation in the present paper, nearly all the masses grow 
with g. However, we see evidence for a state that remains near zero mass as g — > 00. 

The free theory can be solved analytically, and the results are shown in Fig. [H 
Its free-meson spectrum, Fig. Q](a), has K — 1 massless states in each symmetry sector. 
Each such state is made out of two fundamental partons. All states in the corresponding 
glueball sector, Fig. QJb), are massive. The mass scale for all states is set by the CS 
coupling k. To obtain the spectrum, we consider sets of free partons that form mesonic 
color-singlet multi-parton combinations. There are many other combinations of such 
free partons that belong to the non-singlet sector of the Hilbert space, which we can 
omit in the large- N c limit; for g 7^ the only viable states are color-singlets. In other 
words, the free color-singlet combinations are expected to become bound states as soon 
as the coupling is turned on. Mesonic multi-parton color-singlet states are of the form 

fKmi)fl{m 2 )\Q), //(mOa^nO^WIO), //(mOa^nOaLK)//^)^), . . . , 

where the operators /' , p create fundamental partons, while a' create adjoint partons. 
In our theory we have four types of fundamental partons and two types of adjoint 
partons. In the large- N c limit we can have the following four types of mesonic multi- 
parton color-singlet states: 

4(adj.) ii 4|0),ct(adj.) ii 4|0),4(adj.) ii 4|0),4(adj.)^|0), 

where (adj)^ can be any string of adjoint partons. Only the adjoint partons contribute 
to the mass of a state and do so proportional to the CS coupling k. Note also that 
the first pair of fundamental partons above forms a massless combination; there are 
4(K — 1) of these at each value of the resolution. 



S 



To construct the spectrum and find the degeneracies of each mass state, we utilize 
combinatorics and the DLCQ multi-particle formula 



Jmax -* 

1=1 1 

Here 1 < j max < {K — 2) is the number of adjoint partons. Specifically, we calculate the 
compositions Cf = if^Zi) (i- e - ordered partitions) of the integer K into j = (2 + j max ) 
parts, where the factor two counts the two massless fundamental partons in a particular 
state. For example, at K = 5 the compositions that give massive multi-particle states 
are: {Cf,Cf,Cf}. The first composition in the list, Cf, yields (modulo parton type) 
three states with Mf = 5k 2 , two states with M\ = |k 2 , and one state with Mf = |ac 2 . 
Thus, taking into account the several types of partons that can form Cf, the total 
number of states we get for this case is 2 3 (3 + 2 + 1) = 48. The total number of states 
- including the massless ones - as a function of K is thus 

K 

N(K) = ^ k Cf. (3.2) 

k=2 

The dimension N(K) of the Hilbert space of states grows exponentially with K, e.g. 
iV(16) = 28,697,812. 

As an example, consider a four parton state, consisting of two fundamental and two 
adjoint partons, with mass 

M 2 = K, n,m = 1,2,3, .... 

nm 

The two adjoint partons have mass unity (n 2 = 1), while the two fundamentals are 
massless. Thus at K = 7, we have four-parton states starting at M 2 = y. 

The glueball spectrum is evaluated in a similar fashion. However, the free glueball 
multiparticle color-singlet states form closed loops made out of adjoint fermion (fat.) 

and boson (aj.) partons, with a mass easily obtained by (13.11) . Obviously, there are no 
massless states in the glueball sector. The cyclic symmetry of the color trace reduces 
the total number of states that are available, so the free mesons will dominate the 
free energy. Due to supersymmetry it suffices to count only the fermions. Using the 
combinatorics above, we arrive at the number of fermionic states with j partons |24j . 

The function Cf is defined recursively as 

g=l x ' 

Note that C/^^-j-; ^ s zero ^ none of its arguments is an integer. The total 

number of states at a specific K is found by summing over the number of partons j. 
For example, at K — 5 and j = 3, we have a total of 16 states. Eight of these states 
have mass squared M 2 = 10k 2 and the other eight have M 2 = |k 2 . 



9 




20. 



15. 



10. 




GO 



Figure 1: Mass spectra for (a) mesonic bound states and (b) glueball bound states as a 
function of the inverse resolution for 3 < K < 23 when k = 1 and g = (free theory). 
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(a) 



(b) 



Figure 2: Mass spectra for (a) the meson (P-even, T-even) sector and (b) the glueball 
(P-even, T-even) sector as a function of the inverse resolution for 3 < K < 16 when 
k = and g — 1. 
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Figure 3: Mass spectra for (a) the meson (T-even) sector and (b) the glueball (T-even) 
sector as a function of the inverse resolution for 3 < K < 16 when k = 1 and g — 1. We 
note also a small splitting in the masses due to the presence of CS term which breaks 
explicitly the V symmetry. 



3.2 Comparison of meson and glueball spectra 

The generic meson and glueball spectra for nonzero g are shown in Figs. [2] and [H In 
the glueball sector with finite coupling but vanishing CS coupling, there are 2(K — 1) 
massless BPS states. The number of partons in these states grows with the resolution 
K, and there is a mass gap between these massless states and the lowest massive states 
that decreases with increasing resolution. When the CS coupling is not vanishing 
(g 7^ 0, n ^ 0), the BPS massless glueball bound states become approximate BPS 
states [23], with bound-state masses nearly independent of the gauge coupling. The 
masses of the remaining states in this sector grow rapidly with the coupling. 

In earlier work [7] we studied the thermodynamics of this sector with vanishing 
CS coupling. Here, we are considering the mesonic sector of this SYM theory. From 
a previous work [25j I2S] we know that for non-zero coupling there is a mass gap in 
the low-mass sector. The low-mass sector consists of the states that become massless 
bound states of two fundamental partons in the limit that the coupling goes to zero. 
This mass gap decreases as the resolution increases. Of the K — 1 massless states in 
each symmetry sector at vanishing coupling only one remains at finite coupling. 

In the large-iV c limit, the mesonic and glueball sectors decouple. The thermody- 
namics of the theory is generated by the partition function which is the product of the 
partition functions of the two sectors, and the free energy is the sum of the two free en- 
ergies. The glueball bound states are closed loops in color space; their cyclic symmetry 
greatly reduces the number of basis states. Therefore, the number of glueball states 
relative to the meson states at a particular K is very small. One would thus expect 
the mesonic bound states to dominate the thermodynamics. 

There are several ways that the glueball bound states may affect the thermodynam- 
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ics of the full theory. At very low temperature, the thermodynamics will be dominated 
by the very low mass states. At small coupling there are many more light mesonic 
states than glueballs. At strong coupling there are many approximate BPS glueball 
bound states, while only one of the mesonic states remains massless. Thus at strong 
coupling and at temperature high enough to be influenced by the approximate BPS 
states, the thermodynamics will eventually be dominated by the glueball sector. 



4 Density of states and Hagedorn temperatures 

We calculate thermodynamic quantities from the partition function, which we express 
as a sum of Boltzmann factors weighted by the density of states (DoS), p(M 2 , K). The 
discrete spectrum is estimated numerically, and a fit to the data is used to calculate the 
DoS. The discrete spectrum is used to calculate the cumulative distribution function 
(CDF), N(M 2 ,K), which is the number of states with mass squared below M 2 at 
resolution K. The DoS is related to the CDF by 

p(M 2 , K) ee p K (M 2 ) = (4-1) 

with dimensions of L 2 . 

A comment regarding the units of the invariant mass squared eigenvalues M 2 is in 
order. From (12. 16[) it is inferred that the Hamiltonian is of the form 

P- = g 2 A n + K 9 B n + n 2 C n = k 2 { 9 - 2 A n + 9 -B n + C n ) , (4.2) 

where g stands for g^J N c /n. Thus it is a function of a dimensionless ratio, g/n, and 
the dimensionfufl parameter n. The latter sets the mass scale. Here we simply fix the 
value of k to unity, while we numerically investigate the spectra for several values of g. 
So the quantities we calculate are expressed in units where k = 1. 

It is interesting that for g large and k = we have M 2 = g 2 Ai, so the eigenvalues 
scale with g 2 . Therefore, we may determine the strong coupling properties of the theory 
from the solution of the g = 1, k = theory, a numerically much simpler problem. For 
some associated numerical results, see Sec. 15.31 



4.1 Estimating the density of states 

The DoS can be estimated by diagonalizing P~ , computing the CDF from the spectrum, 
and differentiating a smooth fit to the CDF. This is what was done in previous work [7]. 
The size of the matrix representation of P~ increases with K and with the number 
of fields. Eventually, the computational cost becomes too high. To ameliorate the 
situation, we adapted a Lanczos-based algorithm to estimate the CDF directly. 

4 Choosing the CS coupling to set the mass scale is quite natural for the problem at hand, since 
among others we investigate the case where g = and K = 1, where the mass squared eigenvalues are 
proportional to k 2 , sec (|3.1[) . In general we have the freedom of choosing the parameters such that 
they suit the problem. 
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We start by writing the density of states as 

p(M 2 ) = J2dJ(M 2 -M^ (4.3) 

n 

where d n is the degeneracy of the mass eigenvalue M n . The CDF is just 

i-M 2 

Admass 2 < M 2 ) = / dM 2 p(M 2 ). (4.4) 
The density can be written in the form of a trace over e~ lP x+ as follows: 

P^ 2 ) = E dJ(M 2 /2P+ - P~) = -L- jT e^ x+ ' 2P+ d n e-^* + dx+ 

oo 

e iM* X +/2P+ Jre -iP- X + dx+ (45) 



1 



4ttP- 

To approximate the trace, we use an average over a random sample of vectors 
Define a local density for a single vector \s) as 

1 f°° 

Ps(M 2 ) = ^j e iM2 * + ' 2P+ (s\e- ip -* + \s)dx + , (4.6) 
so that the average can be written 

1 S 

p(M 2 )^-J2Ps(M 2 ). (4.7) 

The sample eigenstates \s) can be chosen as random phase vectors [28J, meaning that 
the coefficient of each Fock state in the basis is a random number of modulus one. 

The matrix element (s\e~ lP x+ \s) can be approximated by Lanczos iteration [29]. 
Let D be the square of the norm of |s), and define \u\) — -j^\ s ) as the initial Lanczos 
vector. Then we have 

Ps(M 2 ) = ^f e^ + ! 2P+ {u 1 \e-^* + \u l )dx\ (4.8) 

and (ui\e~ lP x+ \ui) can be approximated by the (1,1) element of the exponentiation 
of the Lanczos tridiagonalization of P~ . Let P~ be this tridiagonal matrix, and solve 
the eigenvalue problem 

M 2 

P7K = i^K- (4-9) 

A diagonal matrix A is related to P~ by the usual similarity transformation P~ = 
UAU' 1 , where Uy = (rf), and Ajj = 5^—^. This means that the (1,1) element is 
given by 

(e- iPa ' x+ ) n = E |«)i| 2 e- iM? ^ +/2P+ . (4.10) 
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M 2 

Figure 4: CDF of the free (g = 0) mesonic T-odd sector at K = 13. Crosses 

(boxes) refer to numerical (analytical) calculation of the bound-state masses. 

The local density is 

p s (M 2 )~-^- / e m2x+ l 2P+ y\{c s n ) l \ 2 e- m ^ x+ / 2P+ dx + 

J n 

- i^T E l«)i| 2 27r5(M 2 /2P + - M 2 J2P+) (4.11) 

n 

^^w sn 6(M 2 -M 2 n ), 

n 

where w sn = D|(c^)i| 2 is the weight of each Lanczos eigenvalue. Note that only the 
extreme Lanczos eigenvalues are good approximations to eigenvalues of the original P~ ; 
however, the other Lanczos eigenvalues and eigenvectors provide a smeared represen- 
tation of the full spectrum. 

The contribution to the cumulative distribution function is 

N S (M 2 )= / dM 2 p(M 2 )~J2 w sn9(M 2 -M 2 n ). (4.12) 

J n 

The full CDF is then approximated by the average 

N(M 2 ) ~ - V N S (M 2 ). (4.13) 

s 

In forming the full CDF, one has to decide how to combine theta functions. This is 
done by using the first sample run as a template for values Mf n at which to evaluate 
N. The contributions of the other samples to N at these values are estimated by linear 
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interpolation in cases where the Lanczos eigenvalues M 2 n are not the same as those in 
the first set. Also, in cases where duplicate eigenvalues are generated by the Lanczos 
iterations, only one is included in the template and the associated weights are added 
together. 

The convergence of the approximation is dependent on the number of Lanczos it- 
erations per sample, as well as the number S of samples. Test runs indicate that the 
recommended value [27] of 20 samples is sufficient. The number of Lanczos iterations 
is kept at 1000 per sample; using only 100 leaves errors on the order of 1-2%. 

A check for the validity of this approach is the comparison between the CDFs for 
the free theory, where the analytic solution is available. Figure H] clearly shows that the 
numerical technique introduced here gives a CDF almost identical to the one obtained 
by the analytical calculation. The very few points that appear to be extraneous have 
no impact on the fitting algorithm we use to calculate the fits to the CDFs. 

4.2 Fits to the spectrum 

In an earlier work related to thermodynamics [7], we split the spectrum into low and 
high-mass regions, separated by the mass gap. The bound-state spectrum for this 
problem has similar characteristics. For instance, for small values of the coupling, 
namely g < 1, the mass gap separates the K — 1 nearly massless color-singlet states 
evolving from the massless states of the free theory from the rest of the spectrum. For 
those values of the coupling these states have M 2 < 1. Thus our density of states, 
Pk{M 2 ), is zero in the mass gap (M 2 ,M 2 ), and the CDF has the following generic 
form: 

f iVi(M 2 ,K), M 2 in <M 2 <M 2 
N(M 2 ,K) = ! const., M 2 < M 2 < M 2 (4.14) 

(N 2 (M 2 ,K), M 2 <M 2 <M^ X . 




2.5 5. 7.5 10. 12.5 15. 2.5 5. 7.5 10. 12.5 15. 0.2 0.4 0.6 0.8 1. 1.2 1.4 

M 2 [X70] M 2 M 2 [x7(T J ] 



(a) (b) (c) 

Figure 5: CDF of the T-even (T + ) sector at K — 14 and g = 0.1. Shown are data 
(dots) and a fit to the data: (a) all states in units of 10 5 states; (b) range of masses 
just above the mass-gap in units of 10 2 states; (c) states below the mass gap. 
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0.2 0.4 0.6 0.8 1. 1.2 1.4 

M 2 

(a) (b) (c) 

Figure 6: Same as Fig. [5] but for the density of states, in units of (a) 10 3 states, (b) 10 2 
states, and (c) 1 state. 

We fit the low-lying mass spectrum of the CDF using the following function 

p(K) 

N 1 (M?,K) = Y,a P M? p , (4.15) 

p=0 

while the logarithm of the CDF for higher masses is fit to the following function 

p(K) 

\n[N 2 (M 2 h ,K)} = (x h + a i yexp[-b 1 M 2 h 5 }J2^ P M 2 h p (4.16) 

p=0 

The other parameters in our fit functions are computed using standard non-linear fit 
algorithms. Typical results are shown in Figs. loTTTOl 




5. 10. 15. 20. 25. 30. 2.5 5. 7.5 10. 12.5 15. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

M 2 [x/0] M 2 M 2 



(a) (b) (c) 

Figure 7: CDF of the T-odd (T _ ) sector at K — 16 and g = 0.5. Shown are data 
(dots) and a fit to the data: (a) all states in units of 10 6 states with point of inflection; 
(b) range of masses just above the mass-gap in units of 10 2 states; (c) states below the 
mass gap. The relatively poor fit near M 2 = 10 in (b) does not have a significant effect 
on the results. 
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(b) 



Figure 8: Same as Fig. [7J but for the DoS, in units of (a) 10 4 states, (b) and (c) 10 
states. 



The spectrum exhibits some structure at relatively small M 2 , as we can see from 
Figs. [5](b) and (c). These states dominate the thermodynamics in the range of temper- 
atures at which we perform the free-energy calculations. Noteworthy is also the point 
of inflection in the CDF plot and the peak in the DoS plot. The data beyond this point 
show the effect of the cutoff imposed by the resolution K. 

^From the plots of the CDF and the DoS, and in particular the figures that depict 
the DoS for g > 0.5, one may predict the result for the free energy T. The main 
contribution to T comes from the one nearly massless state of the spectrum. This 
state exists only in the T-even sector, and consequently this sector dominates the 
thermodynamics at low temperatures. Therefore, we do not expect the fit (e.g., see 
Fig. flOTc)) to give us an accurate result for the free energy, because it essentially leaves 
out the contribution from the nearly massless state. We will return to this point when 
we discuss the numerical results for the free energy in Sec. 15.31 
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M 2 [xiO J ] M 2 M 2 



(a) (b) (c) 



Figure 9: CDF of the T-even sector at K — 16 and g = 4.0. Shown are data (dots) 
and a fit to the data: (a) all states in units of 10 6 states; (b) range of masses just above 
the mass-gap in units of 10 2 states; (c) states below the mass gap. 
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Figure 10: Same as Fig. [91 but for the DoS, in units of (a) 10 3 states, (b) and (c) 1 
state. 



4.3 Hagedorn temperature 

We find that the physical spectrum in the theory grows approximately exponentially 
with the mass of the state 

p H (M 2 ) ~exp(M/T H ), (4.17) 

and therefore has a Hagedorn temperature Th ■ The partition function has the following 
general form 

f M - - f-M\ expfM(^-i)] 
Z^j dMpu (M) exp I — j oc T H T y ^ ■ (4.18) 

The partition function Z obviously diverges as T — > Th, and Th sets the region of 
validity for the calculation of thermodynamic properties. Thus Th serves as an upper 
limit for the temperatures we can use to calculate the thermodynamic functions. 

We will calculate Th by fitting the CDF with an exponential function, and then 
determining its exponent as a function of the resolution and coupling. At fixed coupling 
we extrapolate to infinite resolution and obtain the continuum Hagedorn temperature 
as a function of the coupling. We fit the CDF in the region lying above the mass gap 
and below the point of inflection. The number of states below the mass gap is closely 
related to the number of free massless states and therefore not a factor in the Hagedorn 
domain. The number of states above the point of inflection are significantly reduced 
because of the cutoff imposed by the finite resolution and are therefore not useful in 
determining the Hagedorn temperature. 

In Fig. [TIT a) we show a typical fit to the CDF in this region of M 2 . The data is 
fit well with an exponential. The particular figure deals with the T-odd sector of the 
spectrum at coupling ^ = 0.1 and K = 13. A similar behavior occurs for the T-even 
sector and other values of K. It is clear from our data that the spectral CDF and the 
DoS exhibit a Hagedorn behavior. We plot the logarithm of the CDF versus the bound- 
state mass obtained from our numerical calculations. Then we estimate the range of 
the mass values M where the plot is approximately linear, and fit this region to a linear 
fit of the form aM + T u x . The non-linear part of the distribution, for high values of 
M, is cut off because of the finite resolution K. The extrapolated result for the T-even 
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Figure 11: (a) Logarithm of the CDF versus M. The approximately linear part of this 
logarithmic CDF is fit to f(M) = aM + (b) Extrapolated Hagedorn temperature 
for coupling g = with Tj?(0) « 0.52k. 



sector is very similar to the result of the T-odd sector, so we take the average of the 
two values, for each K. 

In Fig. [TiTb) we show the Hagedorn temperature as a function of the inverse res- 
olution for the massive, free theory (g = 0, n = 1). The data appear to be follow- 
ing a straight line and have been extrapolated to the continuum, where we found 
Tu(g = 0) ~ 0.52 in units where k = 1. We show the Hagedorn temperature in Fig. [12] 
for several cases where (g ^ 0, k — 1) and for resolution K e [11, 16]. As a check we 
note that for the cases g = and g = 0.1 one expects the corresponding extrapolated 
Hagedorn temperatures to be comparable. 

For values of g e {0, 4.0} that are considered here, the upper bound for temperatures 
is set by the Hagedorn temperature of the free theory, T^(0) ~ 0.52k and thus we 
calculate the thermodynamic properties of the theory below this limit. From Fig. [12] 
we glean that Tff(g) grows with the coupling. For larger values of g we may therefore 
access a significantly larger region in T. However, we will leave the discussion of these 
cases for future work. 



5 Finite temperature results in 1 + 1 dimensions 

5.1 The free energy 

We now introduce the basic formulation necessary for our finite temperature calcula- 
tions. Note that our approach here deviates slightly from our earlier work [7J, mainly 
in the way the free energy and the mass squared are normalized. We consider a system 
with constant volume which is in contact with a heat bath of constant temperature. 
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Figure 12: Hagedorn temperature (a) plotted versus 1/K at couplings g = 0.1 (crosses), 
0.5(boxes), 1.0 (triangles), 4.0 (diamonds) in units where k — 1, and (b) extrapolated in 
as function of g with a fit to T^flO = 0.52 + 0.39(7 + 0.054g 2 (dashed line). In (a), 
the dots at 1/K = are the continuum values. For clarity we have included only four 
representative values of g in (a). 

The free energyjf] in units where ks = 1, is given by 

F(T,V) = -ThxZ. (5.1) 

For the large-iV c system at hand, the thermodynamics is described by a canonical 
ensemble of non-interacting glueball and meson-like states. The bound states of the 
theory constitute a supersymmetric two-dimensional free gas. The canonical free energy 
for such a gas in D space-time dimensions is given by 

^^E/^K 1 --^) (6-2) 
^/ = - r E/(|p?rMl + e-*V^) (5.3) 

for bosons and fermions, respectively, where M% in the expression for Th [Th) is the 
invariant bosonic (fermionic) mass spectrum. The integral is performed by expanding 
the logarithm and using the integral representation of the modified Bessel function of 

5 Although our focus is the free energy, T , it is straightforward to use our numerical techniques 
to calculate other thermodynamic functions, such as the internal energy £(T,V) = T 2 (J^§jP) yjthe 
entropy S = {£ - T)/T, and the heat capacity C V (T, V) = (f§) v , cf. [7]. 
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the second kind K„(a;), to find 

OO OO , , , rp\ D/2 



^ M = -2V D .^^) K D/2 (^), (5.4) 



n=l q=l 

/ 71 r m\ D/2 



n=l g=l ^ ^ 7 



(5.5) 



Here Vd-i is the volume in D — 1 space dimensions. The total free energy is obtained 
by adding these two expressions. Because the spectrum is supersymmetric, the sums 
over masses M n traverse the same spectrum, and the total free energy takes the form 

For our calculations we use a rescaled form of Eq. (15.61) . with D = 2 



n=l o=0 



4(K-1)L 2(K-l)^^(2g+l)7r ^ T 

- w^T) h h w^ Kl {—^) • (5 - 7) 

In the last line we introduced the factor dj, which counts degeneracies of mass eigenval- 
ues. This equation is most efficient in the present calculation, because it expresses the 
free energy solely as a function of the numerically evaluated bound-state masses M^. 
We have chosen to normalize by (— 4(K — 1)L) _1 , since 2uq = 4(K — 1) is the total 
number of massless states of the free, massive theory (with g = and n = 1). In prac- 
tice, we can truncate the sum over Bessel functions at q = 10 due to fast convergence. 
Obviously, the sum over states is finite at any finite K. 

The contribution of the massless states to the free energy in D = 2 dimensions can 
be calculated analytically to be 

n b L f°° , p TrTi2 n 



-n,o = -— ^ d Po ^ rF — i =-n,LT- (5. 
•77,0 = r dPozz&TT = -^ LT2 ^ (5-9) 



7T 



for bosons, fermions, and the contribution to the rescaled total, respectively. Thus one 
may separate this contribution from the rest of Eq. (15. 7p . 

Finally, the sum over the states is replaced by an integral over the density of states, 
and Eq. (15. 7h becomes 
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Figure 13: The free energy T for the free theory (g — 0, k = 1) as (a) compared between 
the analytic and numerical results and (b) a function of 1/K. In (a) the temperature 
ranges from 0.015 to 0.5, in units where k = 1, by steps of AT = 0.015k. The solid 
line represents an exact match, the dotted line the actual relation between J- spec t. and 
J 7 fit- In (b) the temperatures are T = 0.1k (crosses), 0.3k (boxes), 0.5k (triangles). 

where p(M 2 ) = pb(M 2 ) = pf(M 2 ) for supersymmetric systems. The T symmetry splits 
the bosonic and fermionic sectors into halves. Thus when calculating the free energy, 
or other thermodynamic properties, we can write 

T = (JS + T f ) r + + {f h + Tj) r - = ftotr+ + FtotT- ■ (5.12) 

5.2 An analytic result: The free energy for the free theory 

Let us start by exploring the free, massive theory [g — 0, k — 1,) which can be solved 
analytically. We will compare the contributions of the meson and glueball sectors to 
the free energy. In particular, using the free-meson sector we can check the validity of 
our approach to replace the sums over discrete spectra with a density of states function 
and check how good our numerical results are compared to an analytic calculation. 

First, we compare the free energy obtained by the means of the analytic method 
outlined above, T spec t., to the free energy extracted from the numerical approach, J- fit, 
using the DoS. The graph presented in Fig. [TBT a) compares analytic and numerical 
results at K = 16 for temperatures 0.015 < T < 0.5, in units where k = 1. We 
deduce from the plot that the agreement is within 1%, which is a typical result. Cutoff 
dependence is very mild: ^F S pect./^fit(K = 13) = 1.012, while ^F spec t./^Ffit(K = 16) = 
1.015. In Fig. [TBT b) we show the free energy of the free, massive theory as a function 
of the inverse resolution at different temperatures. It seems that the free energy T 
converges for low and intermediate temperatures, while it diverges for temperatures 
close to the Hagedorn transition in the continuum limit, as expected. 

We can extract the contributions to the free energy of different parts of the spectrum 
by using Eq. (15.71) . It is interesting to compare the contributions of the two non- 
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interacting sectors of the theory. At temperature T = 0.5k and K = 5, the three- 
parton meson states contribute 1.9 x 10 _2 k 2 to the free energy, while the corresponding 
K = 5, three-parton glueball state contributes only 3.9 x 10™ 4 k 2 . Results at different 
temperatures are listed in Table [TJ The free energy F mesons associated with the free 
meson sector dominates the corresponding F g i ue b a ii- This is a consequence of the fact 
that the mesonic spectrum has 4(K — 1) massless states that contribute |T 2 to F and 
that the meson sector possesses many more states than the glueball sector, especially 
low-mass states which are important in the present calculation. 



T/k 


0.20 


0.25 


0.30 


0.35 


0.40 


0.45 


0.50 


0.55 


F meson j ^ 
F glueball 


0.009 
0.000 


0.016 
0.000 


0.026 
0.001 


0.040 
0.004 


0.059 
0.011 


0.089 
0.026 


0.133 
0.053 


0.199 
0.101 



Table 1: The free energy as a function of the temperature in the meson and the glueball 
sectors in the free theory. 



5.3 Numerical results for nonzero coupling 

We now discuss numerical results for the free energy at finite values of g. These results 
were obtained by applying the DoS method, i.e. replacing the sum over states by an 
integral over the bound-state masses times a density of states computed from a fit 
to the numerically obtained CDF. First, we discuss the temperature dependence and 
then the coupling dependence of the free energy F, for temperatures that lie below the 
zero-coupling Hagedorn temperature Th(q = 0) ~ 0.52k. 

5.3.1 Temperature dependence of the free energy 

The DoS method works also for the interacting theory (g > 0). Using either the discrete 
spectrum approach (sum over the states; (15. 7p ). or the DoS fit to the spectrum yields 
good agreement, at least for weak to medium couplings. The disagreement between 
the two approaches is typically below 1%, as seen in Table [21 This table also presents 
contributions to the free energy from about a thousand states up to M 2 = 9.10283k 2 , 
which belong to the T-even sector at resolution K = 14 and weak coupling g — 0.1. 
We also consider, in Table [3], the free energy at resolution K = 16 for large coupling 
9 = 4.0. 

It is verified from these tables that, at low temperatures, the major contribution to 
the free energy comes from the low-lying states, i.e. the states below the mass gap. In 
the case of Table [3] we have ten such states. As we increase the temperature, more 
states from above the mass gap will contribute significantly. The results are shown in 
Fig. [14] for small coupling, g = 0.1, and in Fig. [15] for large coupling, g = 4.0. For 
K = 16 we find that, at a coupling value of g = 0.1, the free energy is quadratic in T. 
This is expected since for low values of temperature the quasi-massless modes dominate 
and their contribution should be similar to the massless states of the free theory. The 
fourth column of Table [2], which considers contributions from states below the mass 
gap (e.g., M 2 e [0.0001, 0.00113]k 2 ), shows clearly that these states dominate. 
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,07109 


0. 
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1 


,44823 


1 


,44337 


1. 


30232 
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0.3 


1 


,81661 


1 


,81078 


1 


,55613 


0. 


,1359 


0.425 


5 


,33793 


5 


.33113 


3, 


16396 


0. 


2728 


0.45 


6 


,54702 


6 


.54176 


3, 


,55327 


0. 


3059 


0.475 


7 


,98870 


7 


,98597 


3, 


,96515 


0. 


3408 


0.5 


9 


,69247 


9 


.69351 


4. 


39962 


0. 


3776 



Table 2: Free energy as a function of temperature T at K = 14 in the T-even sector 
for weak coupling, g = 0.1: T svect . is obtained by summing over the eigenvalues in the 
interval M 2 e [0.00001, 9. 10283]k 2 ; F fit is obtained by the DoS method described in 
Sec. 14.11 Tf it and JF are the contributions to the latter of the states below the mass 
gap (i.e., M 2 < 0.00113k 2 ) and of a single supersymmetric massless state, respectively. 



At large coupling (g = 4.0) we obtain free energies that are about a hundred times 
smaller than the available weak coupling free energies. This is less than the free energy 
JF that would be contributed from a single massless state! However, it can be justified 
from the spectrum at these temperatures. Recall Figs. M and [IUJ especially those plots 
depicting the spectrum below the mass gap. For the example shown in Table [31 the 
dominant symmetry sector is T-even, the sector which includes the lightest state. We 
denote the contribution of this lightest state by T\. This single state accounts for almost 
all the contribution to the free energy for < T < 0.5k. The latter result is close to the 
contribution JF , that would be made by a single, exactly massless mode, suggesting 
that this very light state may be approximated with a massless state. Finally, at this 
coupling and K = 16, we show in Fig. HBT b) the behavior of the free energy at low T, 
which appears to be quartic. 

5.3.2 Coupling dependence of the free energy 

The behavior of the free energy as a function of the coupling is summarized in Figs. [ToT - 
[TBI For relatively low temperatures (T w 0.1) and for values of g on the order of one 
and above (see Fig. [TBI , the DoS fit misses the most important contribution, which 
is expected from the single lightest state in the T-even sector. For instance, for the 
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Figure 14: Free energy T at weak coupling (g = 0.1) as a function of T at K — 12 
(crosses), 14 (boxes), 16 (triangles) for (a) all temperatures T < T H (dashed vertical 
line) and (b) low temperatures with a quadratic fit F(T) = a w T 2 . 
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Figure 15: Same as Fig. [TU but for strong coupling (g = 4.0). In (b) the low tempera- 
ture behavior is described by a quartic fit F(T) = a s T 4 . 
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[k\ [10" 2 k 2 ] [10~ 2 k 2 ] [10~ 2 k 2 ] [10" 2 k 2 ] 





0.025 
0.05 
0.075 
0.1 






0.0002 
0.0011 
0.0033 






0.0002 
0.0011 
0.0033 






0.0002 
0.0011 
0.0033 





0.0008 
0.0033 
0.0074 
0.0131 



0.225 
0.25 
0.275 
0.3 



0.0384 
0.0505 
0.0644 
0.0803 



0.0383 
0.0503 
0.0641 
0.0796 



0.0383 
0.0503 
0.0639 
0.0793 



0.0663 
0.0818 
0.0990 
0.1178 



0.425 
0.45 
0.475 
0.5 



0.1959 
0.2283 
0.2647 
0.3056 



0.1867 
0.2149 
0.2458 
0.2795 



0.1813 
0.2068 
0.2340 
0.2628 



0.2364 
0.2651 
0.2953 
0.3272 



Table 3: Results for free energy as a function of temperature T at K = 16 and strong 
coupling, g = 4.0. T corresponds to the overall free energy including both symmetry 
sectors. The third column shows the overall contribution of the T-even sector, the 
fourth the contribution from the single nearly massless state, M 2 w 0.0362k 2 , and the 
last column is the contribution that a supersymmetric massless state would make, if it 
were present. 



resolution K = 16 at coupling g = 4.0, the lightest bound state has M 2 = 0.0362k 2 and 
the next available state is at M 2 = 4.86k 2 . Although the fit in Fig.[9](c) seems to capture 
quite well the behavior of the states below the mass gap, states which are expected to 
dominate the thermodynamics at low temperatures, it yields TiT = 0.1) ~ 10 _6 k 2 . 
This is not what we expect from the CDF data. The free energy should be close to 
the contribution of a pair of massless supersymmetric partners, 16( .^_ 1 - ) T 2 . A way to 
improve the calculation of the free energy is to use the discrete spectrum and sum over 
the states instead of approximating this part of the spectrum with a fit function. By 
extracting the states' degeneracies from the CDF data and by utilizing Eq. (]5.7p . we 
obtain TiT = 0.1) ~ 3.33 x 10 _5 k 2 , which matches the expectations much better. In 
fact, the value of T is the contribution of one supersymmetric T-even state. The T-odd 
sector does not contribute significantly to T, since its lightest state (M 2 = 3.651k 2 ) is 
heavily suppressed due to the Bessel factor Ki(M/T) at T = 0.1. 

Although failing here, generally (at relatively weak couplings and small tempera- 
tures) the fit does a good job, mainly because the states below the mass gap are very 
light compared to those for large g, and therefore not suppressed by the Bessel function, 
Kx(M/T), of Eqs. (15. 7p and (15. lip . Therefore, for large values of the coupling, we see 
that as the temperature is gradually being increased, the contribution to the free en- 
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x K=12 
□ K=13 
a K=14 




(b) 



Figure 16: The free energy as a function of the coupling g at temperature T = 0.1k 
(with k — 1) and for resolutions = 12 (crosses), 13 (boxes), 14 (triangles) with two 
different vertical scales. In (b) we see along with the data the contribution to the free 
energy that would be made by a pair of exactly massless superpartners. For fairly large 
values of g and at this temperature, the overall free energy is small compared to the 
contribution of a single pair of massless states. This is expected at this coupling region 
because the masses are very large and are suppressed by the modified Bessel function, 
Ki(x). 



ergy becomes similar to the one coming from an exact massless mode; this contribution 
is included as dotted lines in Figs. [TBI— [151 These results are also in accord with the 
results presented in Tables [2] and [31 

As a check, we have compared results for the massive, strongly coupled theory (g 
large, k — 1) and the massless theory (g = 1, n = 0), where g is the only scale factor. 
We expect the strongly coupled theory to have masses M related to the masses M* of 
the massless theory by M 2 = # 2 M 2 . At g = 4.0, we find M 2 « 4.05 2 M 2 . The free 
energies are related by 

F(T, M\ K) = g 2 F*(T/g, M 2 , K). 

We also calculated the free energy with the DoS method described earlier, and we found 
that it matches quite well the free energy of the theory with g = 4.0 and k = 1. This is 
shown in Fig. |T97 a). Therefore, by solving a numerically less challenging problem, i.e. 
the model with no CS term, we were able to determine the strong coupling behavior of 
the theory with a CS term. 

Having established that g = 4.0 is a relatively strong coupling, and by knowing 
the exact, weak-coupling (g = 0) free energy, let us calculate the strong/weak coupling 
free-energy ratio r s _ w at K = 16, the highest available resolution in our calculations. 
At low temperature, T = 0.1k, we get r s _ w (K = 16) ~ 8.47 x 10~ 3 , and at T = 0.5k we 
obtain r s - w (K = 16) ~ 9.46 x 10 -3 . Results for several temperatures are summarized in 
Fig. [T9T b). A quartic fit to the strong-coupling data reveals that, at resolution K = 16, 
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Figure 17: Same as Fig. [TBI but for T = 0.3. From (b) it is clear that for values of 
g > 1.0 and at this temperature, the contribution to the free energy resembles the one 
from the nearly massless state in the T-even sector. 
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Figure 18: Same as Fig. [161 but for T = 0.5k. From (b) it is clear that more states 
contribute to the free energy at higher temperature, and at relatively large values of g. 
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Figure 19: (a) Comparison of the free energies f{g — 4, k = 1) and .F*(<7* = 1, k* = 0) 
at various temperatures between 0.015k and 0.5k with steps of AT = 0.015k, for 
k = l. The dashed line represents a perfect match, the solid line the best linear fit 
T{g — 4, k — 1) = 0.97.F* (#* = 1,k* = 0). (b) Strong/weak coupling ratio r s - w as a 
function of the temperature at K = 16. 



the ratio is 

r s _ w = -J- « 0.85T 2 . 

This is consistent with the fact that at low temperatures and weak coupling, g, the 
massless states dominate and make a contribution proportional to T 2 to T w . Our CDF 
data suggest that in the continuum limit, however, the lightest (nearly massless) state of 
the strongly coupled theory will become exactly massless, and also yield a contribution 
proportional to T 2 at low temperatures. On the other hand, we know that in the weak 
coupling (g ~ 0) theory for finite K there are exactly 2{K — 1) massless pairs of bound 
states and the ratio r s _ w , for large K, becomes r s _ w (K) = (2(K — l)) -1 . Thus, we 
conclude that in the continuum limit 

K— >oo „ 
Ts—w * "j 

and the discrepancy between the strongly and weakly coupled theories becomes maxi- 
mal. However, we cannot exclude the possibility that we may have more than one pair 
of massless states in the strongly coupled sector of the theory for large values of K, 
namely a number of massless states proportional to K. Although from Fig. [21(a), which 
refers to the strongly coupled system (dual theory), one may try to argue in favor of 
the latter statement that the mass of states at small K seem to follow a trend towards 
the massless limit for relatively large K. However, this is not a definitive result, at least 
from our data, because the highest resolution we have is only up to K = 16. There- 
fore, as far as thermodynamics is concerned in this paper we will just assume that we 
only have one massless pair in the continuum limit of the strongly interacting sector. 
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T r s . w {K = 13) r s - w {K = U) r s _ w (K = 15) r s _ w (K = 16) 





fl0~ 4 l 






fl0~ 4 l 


0.025 


0.12 


0.19 


0.29 


0.50 


0.05 


10.44 


11.32 


12.89 


16.56 


0.075 


42.95 


41.24 


42.71 


49.97 


0.1 


84.74 


76.62 


75.77 


84.72 


1 5 


161 82 


1 38 05 


1 30 57 


1 39 74 


0.2 


217.32 


180.08 


166.48 


174.25 


0.3 


262.02 


208.95 


186.77 


188.81 


0.4 


235.75 


180.94 


154.97 


149.70 


0.45 


207.06 


155.59 


129.71 


121.84 


0.5 


174.62 


128.31 


103.76 


94.57 



Table 4: Data for strong to weak coupling ratio r s - w for K = 13, 14, 15, 16 at various 
temperatures. We show r s _ w (K = 16) as a function of T in Fig. fl9T b). 



Further results for r s _ w for several values of T and K are presented in Fig. [T9Tb) and 
Table HI It seems that r s _ w decreases with K starting at medium temperatures. 

6 Discussion 

We have studied the thermodynamics of M = (1, 1) super Yang-Mills theory in 1+1 
dimensions with fundamentals and a Chern-Simons term that gives mass to the adjoint 
partons. We used SDLCQ to solve the theory in the large-iV c approximation. The 
theory has two classes of bound states: glueballs, which form a closed string in color 
space, and meson-like states, which form open strings. In the large- N c approximation, 
these two sectors do not interact with each other, and make independent contributions 
to the thermodynamics. We previously calculated the contribution of the glueball sector 
but without a CS term. We found that the meson-like sector dominates the glueball 
sector for combinatorial reasons, and, therefore, the results presented here represent 
the full thermodynamics of the theory. Adding a CS term to the theory introduces an 
additional parameter, and thus allows us to inquire about the coupling dependence of 
the theory. 

We have been able to take the calculation up to resolution K = 16, which effectively 
means that we are diagonalizing matrices that are of order 7 x 10 6 by 7 x 10 6 in our 
approximation of the continuum field theory. We introduced a new Lanczos method, 
which is particularly valuable in our calculation of the Hagedorn temperature. 

It is interesting that the spectrum for this theory has a mass gap, which we have 
discussed extensively. The states below the mass gap dominate the low temperature 
behavior of the theory while the states above the mass gap and below the point of 
inflection of the CDF determine the Hagedorn temperature. In fact, the very low 
temperature behavior is dominated by a few massless or nearly massless states in the 
theory. 

The determination of the Hagedorn temperature from the states beyond the mass 
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gap requires a detailed understanding of the SDLCQ spectrum and careful fitting tech- 
niques. We have checked our numerical methods by comparing the solutions of the free, 
massive theory obtained numerically to those extracted analytically. 

We extrapolated the Hagedorn temperature at fixed coupling to the continuum limit. 
The process is repeated at various values of the coupling to determine the coupling 
dependence of the Hagedorn temperature. We find that it increases with the coupling 
from a value of about Tjj = ~k at g = to a value of nearly 3.0k at a coupling of 
9 = 4-0. 

We calculate the free energy of the theory as a function of both the temperature 
and the coupling. As the coupling vanishes, the bound state spectrum can be obtained 
analytically; the analytic results agree with our SDLCQ calculations. The theory has 
4(K — 1) massless fermionic bound states and an equal number of bosonic bound 
states. At low temperature and near-zero coupling, the free energy is simply given by 
the contribution of these massless states, which can be calculated analytically. As the 
temperature increases, the free energy grows quadratically and starts to diverge as the 
temperature approaches the Hagedorn temperature. As we discussed above, this point 
of divergence increases with the coupling. 

At strong coupling and very low temperature, the nearly massless bound states 
dominate the free energy. We find one such fermionic and one such bosonic state. 
These states have very small masses at the highest resolution. Their masses appear 
to decrease with increasing resolution, suggesting that they will become massless in 
the continuum limit. Since the free energy at low temperatures is proportional to the 
number of nearly massless states, the free energy at strong coupling is independent 
of the resolution and therefore has this fixed value in the continuum limit. On the 
other hand, at weak coupling the number of light states grows with the resolution and 
diverges in the continuum limit, as does the free energy. We therefore find that at low 
temperatures the ratio of the free energies at strong and weak coupling goes to zero as 
we approach the continuum. 

A number of interesting extensions of SYM theory with fundamentals and CS term 
exist, for which an SDLCQ mass spectrum can be computed. For instance, one can 
increase the number of dimensions or increase the number of supersymmetries. It would 
be interesting and straightforward to extract the thermodynamic properties of these 
extended theories. 
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